Ultrafast Cascading Theory of Intersystem Crossings in Transition-Metal Complexes 



Jun Chang, A. J. Fedro, and Michel van Veenendaal 
Department, of Physics, Northern Illinois University, De Kalb, Illinois 60115, USA and 
Advanced Photon Source, Argonne National Laboratory, 
9100 South Cass Avenue, Argonne, Illinois 60439, USA 
(Dated: August 16, 2010) 

We investigate the cascade decay mechanism for ultrafast intersystem crossing mediated by the 
spin-orbit coupling in transition- metal complexes. A quantum-mechanical description of the cascad- 
ing process that occurs after photoexcitation is presented. The conditions for ultrafast cascading 
are given, which relate the energy difference between the levels in the cascading process to the 
electron-phonon self energy. These limitations aid in the determination of the cascade path. For 
Fe 2+ spin-crossover complexes, this leads to the conclusion that the ultrafast decay primarily occurs 
in the manifold of antibonding metal-to-ligand charge-transfer states. We also give an interpretation 
for why some intermediate states are bypassed. 



PACS numbers: 82.50.-m 78.47.J- 82.37. Vb 82.53.-k 

I. INTRODUCTION 

Cascade decay is a universal phenomenon associated 
with excited-state relaxation in photophysics and pho- 
tochemistry. In most materials, the excited state that 
is reached after photoexcitation does not directly decay 
to the ground state but follows a complex route of in- 
termediate states with often surprising changes in, e.g., 
spin and lattice parameters. The fastest cascading ef- 
fects are generally associated with dephasing of states 
through the coupling to a continuum, e.g. Fano effects. 1 
After the cascade, the system returns to the ground state 
or to a relatively long-lived metastable state. A highly 
complex example of cascading is photosynthetic water 
oxidation where multiple photoexcitations of a Mn4Ca 
complex bound to amino acid residues in photosystem 
II leads to the production of O2 from water molecules.- 
Obviously, a theoretical quantum-mechanical treatment 
of cascading is of the upmost importance but is also 
highly complex, and even many simpler systems are not 
well understood. A prototypical example of cascading 
occurs in spin-crossover phenomena in transition-metal 
complexes where photoexcitation of a low-spin ground 
state can lead to the creation of high-spin configurations 
on timescales as fast as several hundreds of femtosecond 
(fs). The reverse process has also been observed. Proba- 
bly, the best-studied examples are Fe 2+ complexes with 
a singlet t\ g ground state Q~ A\ in Oh symmetry) and 
a high-spin ^g e g ( 5 ^2) configuration^— The relaxation 
from the metastable high-spin state to the ground state 
is slow with decay times ranging from nanoseconds to 
days. 

The advantage of studying the intersystem crossing in 
Fe compounds is that the cascading clearly involves two 
subsequent ti g ^ — > e g ^ conversions. That this leads to an 
increase in the metal-ligand distance is well known since 
electrons in e g orbitals repel the ligands more strongly 
than those in the t2 g orbitals. However, our understand- 
ing is complicated by the fact that excitations are not 
made into the local dd multiplets, which are generally 



well understood by ab initio techniques in the adiabatic 
limit; 13 ' 14 but into the metal-to-ligand charge-transfer 
states. This makes the exact nature of the cascade path 
difficult to understand due to the competition between 
internal conversion between states of the same spin and 
intersystem crossings between states with different spin. 
Generally, the latter process is considered slower than the 
former. Furthermore, the decay is a nonadiabatic process 
requiring the relaxation of oscillations of vibronic states. 
This intramolecular vibrational energy redistribution is 
considered the fastest process. 

In this paper, we first provide a quantum-mechanical 
cascade decay model for the photoinduced electron state 
in transition-metal complexes. A dissipative Schrodinger 
equation is introduced to include the effects of the in- 
teraction with the surroundings. In the case of Fe spin- 
crossover complexes, we propose a novel and selfconsis- 
tent photon-excited decay path, and find the cascade de- 
cay times in good agreement with experiments on the or- 
der hundreds of femtoseconds from the photoexcited sin- 
glet state to the quintet state. Although our decay times 
are in qualitative agreement with that of phenomenolog- 
ical rate equations, the more detailed understanding of 
the cascading allows a better identification of the states 
involved in the decay path and their time-dependent oc- 
cupation than is possible with rate equations with con- 
stants inferred from experimental data^^ 



II. CASCADE DECAY MODEL 

In a cascading process, several levels at energies Ei are 
involved. The levels couple to the vibrational/phonon 
modes of energy huj of the surrounding ligands. Vari- 
ations in coupling strength \ lead to different metal- 
ligand equilibrium distances for different states. The 
Hamiltonian is given by 

H s = EiJii + fiojo) a + XiUi (a^ + a) , (1) 
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FIG. 1: Left plane: Schematic energy level scheme of Fe 2+ - 
based complexes. The different configurations are given in 
octahedral notation. Manifolds of the metal-to- ligand charge- 
transfer (MLCT) states are indicated by the shadowed areas. 
Right plane: A, e and AR are the energy gap between the 
lowest vibrational levels, the electron-phonon self-energy dif- 
ference, and the change of the equilibrium distance between 
two oscillation states, respectively. 



III. ENVIRONMENTAL DISSIPATION 

The study of the dynamics of dissipative systems is 
notoriously difficult due to the absence of conservation 
of energy and/or particle number in local system and 
the complex interaction with the environment. The den- 
sity matrix method is a standard approach to dissipa- 
tion problems^ However, the density matrices may be- 
come unphysical under the perturbative expansion with 
respect to the system-bath coupling J£ Although Lind- 
blad equations can solve this problem^ a more se- 
rious disadvantage of the density-matrix technique is 
the difficulty in handling large systems. To reduce the 
calculation, a methodology that directly describes the 
nonequilibrium dynamics of the wavefunction is advan- 
tageous. For example, Strunz demonstrated that the dy- 
namics of an open quantum system can be described by a 
non-Markovian stochastic Schrodinger equation. 18 Here, 
we describe a open quantum system with a dissipative 
Schrodinger equation for the system state \ip(t)) which is 
given by 



where rii gives the occupation of state i and is the 
step operator for the vibrational mode. This Hamiltonian 
can be diagonalized with a unitary transformation H s — 
e s H s e~ s , with S = X)i n «v / ft( a ^ — a ) ano - 9i = 
with e l = Af/(7iw), giving, 
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where Hq is the Hamiltonian of the system and D de- 
scribes the effective environmental dissipation. 

When selecting the basis, we write the system vector, 



\m) =5>*(*)w*>' 



(5) 



with eigenstates \ipin) for states i and n excited phonon 
modes. We define the energy difference between the 
states after diagonalization A^- = (E{ — Ei) — (Ej — £j). 
In addition, since a spatial translation of the coordinates 
shifts all the Xi by a constant in H s , only the relative 
change in coupling is of importance. Therefore, it is use- 
ful to define the electron-phonon self-energy difference 
£ij = (A-i — Xj) 2 /(Tilu) between two different states (as 
shown in right panel of Fig. [IJ. An interaction that 
cannot be diagonalized with the electron-phonon term is 
the spin-orbit coupling, which explains the prevalence of 
spin-flips in many decay processes. This causes a cou- 
pling between different states, 



(3) 
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where is the coupling constant and 
particle-conserving transition between states j and i. 
This "local" system is considered part of a larger sys- 
tem such as a molecule in solution or a solid. The lat- 
ter constitute the effective surroundings that can dissi- 
pate energy from the local system. In the following, we 
demonstrate how to incorporate electronic and vibronic 
dissipation. 



with the coefficient Ck in terms of an amplitude ak(t) = 
\ck(t)\ and a phase ipk, or Ck(t) — ak(t)e l(fk ^ . We can 
express the change in the coefficient due to the presence 
of the bath 



dc k {t) 



dt 



da k (t) 



dt 



dt 



.(6) 



In general, the coupling to the environment affects both 
the probability and the phase of the system. The latter 
term gives the change in phase, which causes an embed- 
ding of the local system in its surroundings. Due to the 
complexity of the surroundings, the precise nature of this 
embedding is often very difficult and can usually only be 
taken into account in some effective way. Here, we as- 
sume that the phase of the local system is changed ran- 
domly by the large number of degrees of freedom of the 
surroundings which results in a total phase change close 
to zero according to the law of large numbers. We there- 
fore only consider the changes in the probability by the 
environment. Now let us assume for the moment that we 
are able to determine an expression for the change in am- 



plitude Pfc 



in a particular basis (e.g in the absence 



of certain intersystem couplings) 



dP k 
dt 
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where / is a function of the probabilities Pf.. Below we 
give the explicit expression for the change in probabilities 
related to the damping of phonons. The change in the 
coefficient due to the bath is then given by 



dc k 
dt 



J_ dPk 

2<2fc dt 



1 din Pi. 
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(8) 



This leads to a dissipative term in Eq. Q given by 

_ Ti dhiPk(t) , , . , , . ,„. 



dt 



The dissipation does not necessarily have to be diagonal. 
After deriving the diagonal dissipation in a particular 
basis, a unitary transformation to a different (more suit- 
able) basis can be made. After deriving an expression 
for the dissipation in a particular basis set, we can solve 
the problem in the presence of intersystem couplings and 
dissipation. 

In cascade decay model, the eigenvectors \ipi n ) of H s 
are selected as the basis. The vibrational cooling by 
the bath can be taken into account by the dissipative 
Schrodinger equation in Eq. (j4]) with Hq substituted by 
H a + Hi and Hj = e s Hje~ s . We still need the detailed 
formulae for Eq. ([7]). For state i with n excited phonon 
modes, on the one hand, the vibrational coupling to the 
surroundings relaxes a state with n phonons to a n — 1 
phonon state by the emission of phonons. On the other 
hand, the probability of the n-phonon state increases due 
to the decay of the state with n + 1 phonons i 19 ' 20 This 
gives a change in the probability of the n-phonon state 



dP in (t) 
dt 



= -2nTP m (t) + 2( ? i + l)rp,„ +1 , 



(10) 



where P m (t) = \(i/) in \ i/}(t))\ 2 , T = npV 2 /h is the re- 
laxation constant, where p is the effective bath phonon 
density of states and V is the interaction between the 
local system and the bath. Due to the complication of 
V, we take T as a parameter. Note that the decay time 
from state i to j is not directly related to T. Our numeri- 
cal calculations are not sensitive to the change in (2r) _1 
from 20 fs to 60 fs. In Figs. g]and[3J we take (2F)" 1 = 30 
fs. 



IV. IRON SPIN-CROSSOVER COMPLEXES 

Let us specifically consider iron-based complexes, such 
as Fe 2+ (bpy)3. A typical energy level scheme is shown 
in Fig. [TJ The Fe d 6 multiplet levels have been ex- 
tensively studied. However, since these excitations are 
dipole forbidden, the system is often excited into so- 
called metal-to-ligand charge-transfer (MLCT) states. 
Spin crossovers are generally strongly covalent, where the 
metal states are given by a\d 6 }+J2 i(Pi\d 7 L.i)+Ji\d 5 Li)) + 
• • • , where and are an electron and a hole, re- 
spectively, on the ligand labeled i. The MLCT states 




FIG. 2: (color online) The time t 90 % for the third state 
reaching its 90% population as a function of the energy dif- 
ference between state 2 and 3 with A12 = £12 = 0.2 eV 
and £23 = 0.4 eV. The solid (red) curve is calculated us- 
ing the dissipative Schrodinger equation, the dashed (green) 
curve is based on rate equations with t 90 % = — In 0.1 /70 with 
7o = 7172/(71 +72), 71,2 are defined in Eq. 1121 For the values 
of A23 indicated by the squares, the time dependence is given 
in Fig. El 



are the antibonding states related to the metal states 
(e.g. 1 MLCT- 1 Ai, 3 MLCT- 3 Ti )2 and 5 MLCT- 5 T 2 ). 
Since the photoexcitation conserves spin, the initial pho- 
toexcited state starting from a 1 A\ ground state is the 
1 MLCT state (with predominantly t^L 1 character). A 
major problem in the identification of the states in the 
cascading process is that optical selection rules make 
it impossible to directly probe the intermediate states 
spectroscopically, and several cascading paths to the 
metastable high-spin 5 T2 configuration have been pro- 
posed. McCusker et al& suggested that the 1 MLCT state 
directly decays to the high-spin 5 T2 quintet state requir- 
ing a spin change of AS = 2. Gawelda et al& proposed 
that the 1 MLCT state relaxes to the 3 MLCT within 30 
fs followed by a departure from this state within ^120 fs. 
Subsequently, the quintet state is reached via interme- 
diate ligand-field multiplets. Recently, Bressler et alr^ 
claimed that the 3 MLCT (t^e^L 1 ) state directly relaxes 
to the 5 T 2 state bypassing the ligand field triplet state. 
It is surprising that several states are bypassed, such as 
the triplet states 3 Ti j 2, which are located between the 
3 MLCT and 5 T 2 , and that there is no 3 MLCT to 3 Ti :2 
internal conversion, which is expected to occur on a faster 
timescale than the spin-crossover process. In order to un- 
derstand the cascade process in Fe-compounds, we first 
need to establish an appropriate range of parameters for 
Fe spin-crossover compounds. The strength of the inter- 
action between the iron and its surrounding ligands can 
be obtained from ab initio calculations^ The change in 
energy for different configurations is close to parabolic 
for an adiabatic change in the Fe-ligand distance. From 
the change in equilibrium distance, we can obtain a dif- 
ference in electron-phonon coupling of approximately 
£12 = 0.2 eV between the singlet and triplet state, and 
£23 = 0.4 eV between the triplet and quintet state. A 
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FIG. 3: (color online) The probability of finding a particular 
state as a function of time for different values of A23. The 
green dashed-dotted, the blue dashed, and the red solid curves 
give the probability of finding states 1, 2, and 3, respectively: 
a) A 23 = 0.4 eV; (b) A 23 = 0.09 eV; (c) A 23 = 0.69 eV. Other 
parameters are the same as those in Fig. [2] 



typical energy for the Fe-ligand stretching mode is Tvuj — 
30 meViSi Next, we need to determine the energies of 
the states involved in the cascading process. Since elec- 
tronic transitions do not directly change the metal-ligand 
distance, an estimate of 2.6 eV can be obtained for the 
energy difference between the lowest vibrational levels of 
the 1 MCLT and 1 A\ states from the pump laser wave- 
length of 400 nm and the fluorescence of 600 nm from 
the lowest 1 MCLT state back to the 1 A\ state. Another 
weaker emission shows a shift to 660 nm corresponding to 
a change in energy of A 12 = 0.2 eV4 This small energy 
difference implies that state 2 is a 3 MLCT stated see Fig. 
[TJ For the spin-orbit coupling, we take the atomic value 
for Fe, V12 = V23 = 0.05 eV. This leaves us with the 
question of the nature of the third state and how a com- 
plete relaxation to the lattice parameters of the quintet 
state can be observed within 300 fs. We can obtain more 
information on the energy of the level in the cascading 
process by solving the cascading equations given above 
for cascading from level 1 to level 3 via an intermediate 
level 2. We define i state occupation Pi{t) = J2 n Pin{t)- 
Figure [2] shows the results for tgo%, i.e., the time when 
the third level reaches its 90% population, as a func- 
tion of A23. The value of t 90 % shows a broad minimum 
around A23 = £23 — 0.4 eV. The minimum calculated 
time to reach a 90% occupation of the third level in the 
cascading process is about 200 fs. This is in good agree- 
ment with the results by Bressler et al^ that show that 
the metal-ligand distance reaches the high-spin value on 
the same timescale. Let us first note that £90% increases 
rapidly when A23 is less than 0.3e23 = 0.12 eV or larger 
than 1.7623 — 0.68 eV. This implies that the total en- 
ergy bridged in the fast cascading process is at most 
A12 + A23 = 0.32-0.88 eV. This is insufficient to over- 
come the more than 1.8 eV energy gap between the initial 
photoexcited state and the 5 T2 high-spin state. Another 



possibility is that an internal conversion occurs between 
states 2 and 3, for example, state 3 is a metal-centered 
3 ?i i 2 triplet state, which are about 1 eV below state 2. 
For an internal conversion, the electron number of e g or- 
bitals is unchanged, and £23 has to be less than £12 which 
does involve a t2 g ^ — > e g \. conversion. Using a reduced 
electron-phonon self-energy of £23 of 0.2 eV while keep- 
ing V23 = 0.05 eV, comparable results are obtained for 
the cascade time as a function of A 2 3/£23 as in Fig. [2] 
(not shown). In order for the decay to occur in less than 
a picosecond, the energy gap between states 2 and 3 has 
to be less than ~2e23 = 0.4 eV. Therefore, even though 
the energy gap for an internal conversion is smaller, so 
is the maximum energy difference for which we can have 
ultrafast decay, and can also rule out the metal-centered 
triplet states as the third state in the cascading process. 
This is in agreement with experiment, since the metal- 
ligand separation increases to 0.2 Ai . showing that the 
high-spin state is reached. Furthermore, many internal 
conversion processes are also bypassed. Since the ligands 
7r* orbitals are almost orthogonal to the metal e g or- 
bitals, the electrons at the ligand atoms can only return 
to the ti g orbital. Therefore, this prohibits, for example, 
decay from the 1 MLCT (t^L 1 ) to the ^ state (& g e\). 
The only remaining possibility is that state 3 is a 5 MLCT 
(i^e^L 1 ) state. Furthermore, in the broad 5 MLCT man- 
ifold, it is easier to find a level with A23 ~ £23- This 
leads to the interesting conclusion that the entire ultra- 
fast decay process occurs primarily in the antibonding 
MLCT states, a possibility not considered before. Sub- 
sequently, the 5 MLCT state can relax more slowly to the 
5 T2 state through internal conversion, which is difficult to 
observe in EXAFS experiments^ due to the comparable 
lattice parameters for the quintet states. This relaxation 
is slower since the spin-orbit coupling does not couple 
the two quintet states. Another reason that the ultrafast 
decay occurs in the MLCT states is that the broad man- 
ifold of MLCT states makes it easier to find states for 
which Ay ss £y, whereas for the ligand field multiplets, 
this condition is only satisfied accidentally. Therefore, 
for the first step in the cascade decay, it is not a co- 
incidence that the strongest increase in intensity in the 
fluorescence occurs at A12 « 0.2 eV^ close to the value of 
£12 although the energy of 3 MLCT manifold spans over 
1.5 eV. 

A qualitative understanding of the time dependence of 
the decay can be obtained starting from classical phe- 
nomenological rate equations^^ 



dPjjt) 
dt 



7 i _iP i _ 1 (t)-7 i P,(t), 



(11) 



where Pi(t) is the probability of ith level with i = 1...N. 
The rate constants % can be calculated using Fermi's 
golden rule ) 22 i 23 



7l = 2nF n V l ? i+1 /h 2 iV, 



(12) 



where V^j+i is the interaction between levels and 
F n = e~ 3 g n /n\ is the Franck-Condon factor with n w 
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Aj^+i/Tkd; g = Ei^+i/Tiuj is the Huang- Rhys factor. 
Solving the above equations, we can approximately write 
the occupation of the final level in the cascading 

P N {t) « 1 - e-*>*, (13) 

with I/70 = X^il^ 1 (l/7»)- Figure [2] shows that, al- 
though rate equations give a qualitative description of 
igo% as a function of A23/623, quantitative differences 
are present. In the rate equation, energy conservation 
in Fermi's Golden Rule restricts the relaxation to states 
of equal energy, whereas in the quantum mechanical cal- 
culation, the state couples to the whole Franck-Condon 
continuum. This leads to a smaller minimum value for 
tgo% and also to a significantly broader minimum for 
the quantum-mechanical cascading process. The broad 
width of the minimum somewhat loosens the condition 
Ay rs £23, which further explains the prevalence of 
ultrafast-decay processes. 

For A23 = 0.4 eV, we study the time-dependent oc- 
cupations of the three states involved in the cascading 
process, see Fig. |^a). We find that state 1 decays with 
a relaxation time about 20 fs. While the probability of 
state 2 increases quickly in the first 20 fs, it then begins 
to decrease due to relaxation into state 3. At 120 fs, it 
has lost most of its population, which agrees well with 
the experimentally observed departure from this state 
within ~120 fs£ State 3 reaches almost 100% within 300 
fs. All these timescales agree well with experiments.—^ 
For A 23 = 0.09 eV and 0.69 eV, see Fig. fflb) and (c), 
the decay of state 1 is almost the same as for A23 = 0.4 
eV, since the decay from 1 — > 2 has been left unchanged, 
but the decay from state 2 has slowed down dramatically. 



in transition-metal complexes. Ultrafast cascading oc- 
curs when the energy difference between the levels is 
comparable to the self energy e. Since the latter is on 
the order of several tenths of an electronvolt, restrictions 
are imposed on the energy gap that can be bridged in the 
ultrafast cascading process. On the other hand, the man- 
ifold of the metal-to-ligand charge-transfer states make 
the gaps between energy levels many and various. As a 
result, the ultrafast cascading in Fe spin-crossover com- 
plexes after excitation with visible light occurs primarily 
in these manifolds. We propose a novel and selfconsis- 
tent photon-excited decay path for [Fe n (bpy)3] 2+ , which 
is then 1 MLCT^ 3 MLCT^ 5 MLCT-> 5 T 2 Va x . Good 
agreement is found between the calculated and exper- 
imentally observed decay times. Our quantum model 
results qualitatively agree with the calculation of decay 
times based on the phenomenological rate equations com- 
bining with Franck-Condo factor. We further give an 
explanation why some intermediate states are bypassed. 
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V. CONCLUSIONS 



In conclusion, we have provided a quantum-mechanical 
model for the cascade decay mechanism of spin crossover 
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